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O ■ We address the issue of whether a region of disordered chiral condensate (DCC), in 

Qv . which the chiral condensate has components along the pion directions, can form. We 

f-| ' consider a system going through the chiral phase transition either via a quench, or via 

^^. relaxation of the high temperature phase to the low temperature one within a given 

^H. time scale (of order ~ Ifm/c). We use a density matrix based formalism that takes 

,^ . both thermal and quantum fluctuations into account non-perturbatively to argue that 

if the 0(4) linear sigma model is the correct way to model the situation in QCD, then 

^ ■ it is very unlikely at least in the Hartree approximation, that a large (> 10 fm) DCC 

j^ ' region will form. Typical sizes of such regions are ~ 1 — 2 f m and the density of pions 

in such regions is at most of order ~ 0.2/fm . We end with some speculations on how 

large DCC regions may be formed. 



1 Introduction 

The proposition has been put forth recently that regions of misahgned chiral vacuum, or 
disordered chiral condensates (DCC) might form in either ultra-high energy or heavy nuclei 
collisions [|1], 0, ^ Q. If so, then this would be a striking probe of the QCD phase transition. 
It might also help explain 0, Q the so-called Centauro and anti-Centauro events observed in 
high-energy cosmic ray experiments^. 

How can we tell, from a theoretical standpoint, whether or not we should expect a DCC 
to form? Clearly, investigating QCD directly is out of the question for now; the technology 
required to compute the evolution of the relevant order parameters directly from QCD is 
still lacking, though we can use lattice calculations for hints about some aspects of the QCD 
phase transition. What we need then is a model that encodes the relevant aspects of QCD 
in a faithful manner, yet is easier to calculate with than QCD itself. 

Wilczek and RajagopalQ have argued that the 0(4) linear a-model is such a model. It 
lies within the same static universality class as QCD with two massless quarks, which is a 
fair approximation to the world at temperatures and energies below Aqcd- Thus work done 
on the n = 4 Heisenberg ferromagnet can be used to understand various static quantities 
arising at the chiral phase transition. 

One conclusion from reference(§) was that, if as the critical temperature for the chiral 
phase transition was approached from above the system remained in thermal equilibrium, 
then it was very unlikely that a large DCC region, with its concomitant biased pion emission, 
would form. 

The point was that the correlation length ^ = m~^ did not get large compared to the 
T~^. A more quantitative criterion involving the comparison of the energy in a correlation 
volume just below Tc with the T = pion mass (so as to find the number of pions in a 
correlation volume) supports the conclusion that as long as the system can equilibrate, no 
large regions of DCC will form. 

The only option left, if we want to form a DCC, is to insure that the system is far out of 
equilibrium. This can be achieved by quenching the system (although Gavin and M tiller |I0| 
claim that annealing might also work). What this means is the following. Start with the 
system in equilibrium at a temperature above T^. Then suddenly drop the temperature to 
zero. If the rate at which the temperature drops is much faster than the rate at which the 
system can adapt to this change, then the state of the system after the quench is such that 
it is still in the thermal state at the initial temperature. However, the dynamics governing 
the evolution of that initial state is now driven by the T = Hamiltonian. The system will 



then have to relax from the initial state, which is not the ground state of the Hamiltonian 
to the zero temperature ground state. During this time, it is expected that regions in which 
the order parameter is correlated will grow. We can then hope that the correlation regions 
will grow to be large enough to contain a large number of pions inside them. 

The possibility that the chiral phase transition might occur following a quench in heavy- 
ion collisions was explored by Wilczek and Rajagopal in reference(0). They argue there 
that long wavelength fluctuations in the pion fields can develop after the quench occurs. 
Modes with wavenumbers k smaller than some critical wavenumber fccrit will be unstable 
and regions in which the pion field is correlated will grow in spatial extent for a period of 
time. The essence of this mechanism is that the pions are the would-be Goldstone particles 
of spontaneous chiral symmetry breaking. In the absence of quark masses, the pions would 
be exactly massless when the vf and a fields are in their ground state. However, during the 
quench, the a field is displaced from its zero temperature minimum so that the required 
cancellation between the negative bare mass^ term in the Lagrangian and the mass^ induced 
through the pion interactions with the a condensate does not occur. This then allows 
some of the pion momentum modes to propagate as if they had a negative mass^, thus 
causing exponential growth in these modes. In Ref.(0), the classical sigma model was 
simulated. The correlation functions were taken as spatial averages and the expectations for 
the growth of various pion field momentum modes were borne out. More recently however. 



Gavin, Gocksch and Pisarski[TT[] have concluded that the strongly coupled linear sigma 
models does not produce large correlated domains of pions. These authors also performed a 
numerical simulation of the classical equations of motion. However, there are already hints 



from previous work[0, |13| that both quantum and thermal effects, may be important in 
determining the growth of correlation regions in a field theory. It is to this calculation that 
we devote the rest of this paper. In the next section, we will develop the formalism necessary 
to take the quantum and thermal fluctuations of the fields into account. Having done that, 
we turn to the actual numerical solution of the equations we find, and use these solutions to 
calculate the equal time two-point correlation function for the pions. It will be clear after 
doing this that in the case in which the system is either quenched or relaxed from a high 
temperature phase whose temperature is larger than the critical temperature, the correlation 
regions are never large enough for the correlations in pions to be observed. We strengthen 
this conclusion by computing the number of pions in the correlation volume. We then end 
with some speculations concerning possible ways in which a large DCC region might form. 



2 The 0(4) cr-Model Out of Equilibrium 



Our strategy is as follows. We will use the techniques developed by us previously jT^ and 



use the functional Schroedinger representation, in which the time evolution of the system is 
represented by the time evolution of its density matrix. 

The next step is to evolve the density matrix in time from this initial state via the 
quantum Liouville equation: 

zn^ = [H,p{t)], (1) 

where H is the Hamiltonian of the system after the quench. Using this density matrix, we 
can, at least in principle, evaluate the equal time correlation function for the pion fields, and 
observe its growth with time. 

Let us now implement this procedure. We start with the sigma model Lagrangian density 



£ = -9^<l>-9^$-y((T,7f) (2) 

V{(T,n) = -m2(t)$.$ + A($-$)2-/ia (3) 

where $ is an 0{N + 1) vector, $ = (cr, vf) and vf represents the A^ pions. 

The linear sigma model is a low energy effective theory for an SUi^{2) x SU^{2) (up 
and down quarks) strongly interacting theory. It may be obtained as a Landau- Ginzburg 



effective theory from a Nambu-Jona-Lasinio model|15[. In fact, Bedaque and Das[0 have 
studied a quench starting from an SUi,{2) x S'[/r(2) Nambu-Jona-Lasinio model. 

We have parametrized the dynamics of the cooling down process in terms of a time 
dependent mass term. We can use this to describe the phenomenology of either a sudden 
quench where the mass^ changes sign instanteneously or that of a relaxational process in 
which the mass^ changes sign on a time scale determined by the dynamics. In a heavy ion 
collision, we expect this relaxation time scale to be of the order of r ~ 0.5 — Ifm/c. The term 
ha accounts for the explicit breaking of chiral symmetry due to the (small) quark masses. 
We leave A^ arbitrary for now, though at the end we will take A^ = 3. 

Our first order of business is to identify the correct order parameter for the phase tran- 
sition and then to obtain its equation of motion. Let us define the fluctuation field operator 
X(x,t) as 

cr = ^it) + xix,t), (4) 



with 0(t) a c-number field defined by: 

Here Q is the spatial volume we enclose the system in. The fluctuation field xi^^ t) is defined 
so that (i) {x{x,t)) = 0, and (ii) xi^^t) = ~4'if)- Making use of the Liouville equation for 
the density matrix, we arrive at the following equations: 

m = P(^) = ^l d'x{U^{x)) (6) 

where Ilcr{x) is the canonical momentum conjugate to (j{x). 

The derivative of the potential in the equation for 7r(t) is to be evaluated at cr = 0(t) + 
x{x,t). These equations can be combined into a single one describing the evolution of the 
order parameter 0(t): 

0(t) + yd'x (^^^) \.=HtHxi^^t)) = 0. (8) 

To proceed further we have to determine the density matrix. Since the Liouville equation 
is first order in time we need only specify p{t = 0). At this stage we could proceed to a 
perturbative description of the dynamics (in a loop expansion) . 

However, as we learned previously in a similar situation |[1^, |13|, the non-equilibrium 
dynamics of the phase transition cannot be studied within perturbation theory. 

Furthermore, since the quartic coupling of the linear sigma model A must be large (A ~ 
4 — 5 so as to reproduce the value of /tt ~ 95 Mev with a "sigma mass" ~ 600 Mev), the 
linear sigma model is a strongly coupled theory and any type of perturbative expansion will 



clearly be unreliable. Thus, following our previous work[12, H and the work of Rajagopal 



and Wilczek0 and Pisarski|TT|] we invoke a Hartree approximation. 

In the presence of a vacuum expectation value, the Hartree factorization is somewhat 
subtle. We will make a series of assumptions that we feel are quite reasonable but which, 
of course, may fail to hold under some circumstances and for which we do not have an a 
priori justification. These are the following: i) no cross correlations between the pions and 
the sigma field, and ii) that the two point correlation functions of the pions are diagonal 



in isospin space, where by isospin we now refer to the unbroken 0{N) {N = 3) symmetry 
under which the pions transform as a triplet. These assumptions lead to the following Hartree 
factorization of the non-linear terms in the Hamiltonian: 



X^ ~^ 6(x^) + constant (9) 

X' - 3(x')x (10) 

(tt-tt)^ -^ (2 + —)(7f 2) ^2 ^constant (11) 

^V - ^'{X') + {^')X' (12) 

v?'x ^ (vf')x, (13) 

where by "constant" we mean the operator independent expectation values of the composite 
operators. Although these will be present as operator independent terms in the Hamiltonian, 
they are c-number terms and will not enter in the time evolution of the density matrix. 

It can be checked that when = one obtains the 0{N + 1) invariant Hartree factor- 
ization. 

In this approximation the resulting Hamiltonian is quadratic, with a linear term in x'- 



^„,)./.».{S,|,(^,(^,,V-W.^X=.^.| (14) 

Here H^, Htt are the canonical momenta conjugate to xi^)^ t^{x) respectively and V^ is 
recognized as the derivative of the Hartree "effective potential" [0, |T^ with respect to (f) (it 



is the derivative of the non-gradient terms of the effective action fl^, p^ , |T9[)- 

In the absence of an explicit symmetry breaking term, the Goldstone theorem requires 

the existence of massless pions, M.t^ = whenever V^ for (p ^ 0. However, this is not the 

case within our approximation scheme as it stands. 

This situation can be easily remedied, however, by noting that the Hartree approximation 

becomes exact in the large N limit. In this limit, (vf^) ^ 0{N), {x^) ^ C(l), 0^ ~ 0{N). 

Thus we will approximate further by neglecting the 0{1/N) terms in the formal large A^ 

limit. This further truncation ensures that the Ward identities are satisfied. We now obtain 

V\t) = 0(t) [m2(t)+4A0'(t)+4A(7f2)(i)] -/i (15) 

MUt) = m=^(t)+4A</)2(t)+4A(7f2)(t) (16) 

Mlit) = m\t) + 12\<p\t)+A\{n^)it). (17) 
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The Hamiltonian is now quadratic with time dependent self-consistent masses and Gold- 
stone's Ward identities are satisfied. 

Since the evolution Hamiltonian is quadratic in this approximation, we propose in the 
Hartree approximation a gaussian density matrix in terms of the Hartree-Fock states. As a 
consequence of our assumption of no cross correlation between x ^-i^d tt , the density matrix 
factorizes as 

p{t) = p^{t) ® p^{t) (18) 

In the Schroedinger representation the density matrix is most easily written down by 
making use of spatial translational invariance to decompose the fiuctuation fields x{^) t) ^-i^d 
7r(x, t) into spatial Fourier modes: 



Xix,t) 



1 



Y,Xk{t)Gw{-ik- 



X) 



TT[X) 



^7f^exp(-i/c 



x) 



(19) 
(20) 



where we recall that xi^^ t) = (^{x) — 0(t). We can now use these Fourier modes as the basis 
in which to write the density matrices for the sigma and the pions. We will use the following 
Gaussian ansatze: 



Px[x,X,'t] = n^x,fe(^)exp 



^"''^'^■x.(t)x-.(t) + %^x.(t)x-.(t)+ 



2h 



2h 



Bx,kit) 

h 



Xk{t)x-k{t) 



+ j^PxAt){x^k{t)-x-k{m, 



(21) 



p^ [tt , TT , t] = Y[ ■^^,k (t) exp 



■TTfc ■ 7r_fe H — TTfe ■ 7r„fc + 



2h 



2h 



Bjrjtit) 



-TXk ■ TX-k 



(22) 



The assumption of isospin invariance implies that the kernels A^ fc, 5^ ^ transform as isospin 
singlets, since these kernels give the two point correlation functions. Furthermore, hermiticity 
of the density matrix requires that the mixing kernel B be real. The lack of a linear term in 
the pion density matrix will become clear below. 

The Liouville equation is most conveniently solved in the Schroedinger representation, in 
which 



7 



n^(f) 



-ih-— 

ox 



Uiix) 






ih 



dp{t) 
dt 



H[U^, n^, X, vr; t] - H[I{^, H^, x, vf; t] p(t) 



(23) 



Comparing the terms quadratic, linear and independent of the fields (x ; tt), we obtain 
the following set of differential equations for the coefficients and the expectation value: 



.X 



x,fc 



K 



x,k 



lA 
iB 



x,k 



X,k 



^l,kit) 



-(^x,fc-^x,A 



_K,k - B^k 


2 


it)_ 


B^,k \A^,k - 


^xj 




e + Mlit) 







4> + m^{t)4> + 4A0^ + 4A0(7f2(f , t))-h = 0. 



(24) 

(25) 
(26) 
(27) 

(28) 



.K 



TT,k 






;(^x,A 



^x,fc) 



<,-5^,,-<,(t) 

-Djr.fc ( ^7r,fc ~ ^TT,k 
2 I A ^2 



fc^+^^(t). 



(29) 

(30) 
(31) 
(32) 



The lack of a linear term in ( P^D is a consequence of a lack of a linear term in tt in 
the Hartree Hamiltonian, as the symmetry has been specified to be broken along the sigma 
direction. 

To completely solve for the time evolution, we must specify the initial conditions. We 
will assume that at an initial time (t = 0) the system is in local thermodynamic equilibrium 
at an initial temperature T, which we take to be higher than the critical temperature, T^ ~ 
200 MeV , where we use the phenomenological couplings and masses to obtain T^. 

This assumption thus describes the situation in a high energy collision in which the 
central rapidity region is at a temperature larger than critical, and thus in the symmetric 
phase, and such that the phase transition occurs via the rapid cooling that occurs when the 
region in the high temperature phase expands along the beam axis. 



The assumption of local tliermodynainic equilibrium for the Hartree-Fock states deter- 
mines the initial values of the kernels and the expectation value of the sigma field and its 
canonical momentum: 

A^^k{t = 0) = u^^kiO) coth[(3huj^^k{0)] (33) 

A.,fc(t = 0) = u^,k{0) coth[phiu^^kiO)] (35) 



4>{t = 0) = 00 ; 0(t = 0) = 0, (37) 

with f3 = l/ksT. We have (arbitrarily) assumed that the expectation value of the canonical 
momentum conjugate to the sigma field is zero in the initial equilibrium ensemble. These 
initial conditions dictate the following ansatze for the real and imaginary parts of the kernels 
Ajr fc(t), A^^kit) in terms of complex functions A-jr^kit) = An-Tr^kit) + iAi-^-jr^kit) and A^^kit) = 
AR;^,k{t) + iAi;y,,k{t)M- 



AR.^,kit) = AR.^,k{t) coth[phu^^kiO)] (38) 

B.At) = - • ufr'^'Li (39) 

Ai.,^,kit) = Ai.,^,kit) (40) 

The differential equation for the complex function A can be cast in a more familiar form by 
a change of variables 

A.At) = -^^r^y (41) 

with \l/^ fc obeying the following Schrodinger-like differential equation, and boundary condi- 
tions 

rf2 



df" 



+ <kit) 



^.,fc(t) = (42) 

1 



^.,fc(t = 0) = ; *^,fc(t = 0)=z^cu^,fc(0). (43) 

yuJn,kiO) 

Since in this approximation the dynamics for the pions and sigma fields decouple, we will 
only concentrate on the solution for the pion fields; the effective time dependent frequen- 
cies for the sigma fields are completely determined by the evolution of the pion correlation 
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functions. In terms of these functions we finally find 



Nh 



(7ffc(t) ■ ^_fe(t)) = ^ I ^^,fc(t) 1^ COth 



2kBT 



(44) 



In terms of this two-point correlation function and recognizing that the \I'7r,fc(^) only depends 
on k"^, we obtain the following important correlations : 



(7r(f,t)-7f(0,t)) 



4^ 
Nh 
4^ 



dkk^ I ^^,A:(^) 

sin(A;x 



coth 



2kBT 



dkk- 



X 



^^,fc(t) P coth 



2kBT 



(45) 

(46) 
(47) 



The presence of the temperature dependent function in the above expressions encodes the 
finite temperature correlations of the initial state. The set of equations ( [28|J^ ) with the above 
boundary conditions, completely determine the non-equilibrium dynamics in the Hartree- 
Fock approximation. We will provide a numerical analysis of these equations in the next 
section. 



2.1 Pion production 

In the Hartree approximation, the Hamiltonian is quadratic, and the fields can be expanded 
in terms of creation and annihilation of Hartree- Fock states 



Mt) = -^ (4^l,fc(t) + aLfc^.,fc(t) 



(48) 



The creation a|, and annihilation a^ operators are independent of time in the Heisenberg 
picture and the mode functions '^TT,k{t) are the solutions to the Hartree equations ( ^2]) which 
are the Heisenberg equations of motion in this approximation. The boundary conditions ( ^3|) 
correspond to positive frequency particles for t < 0. The creation and annihilation operators 



may be written in terms of the Heisenberg fields (|48|) and their canonical momenta. Passing 
on to the Schroedinger picture at time t = 0, we can relate the Schroedinger picture operators 
at time t to those at time t = via a Bogoliubov transformation: 



Skit) = J-+,fc(t)afc(0) + J^-,fe(t)al,(0) 



(49) 



with 
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2 1 I ^^,fcW 



4|^.,fc(0) P 



\ ^ I ^^,fcW 



<,(0) I vl/,,,,(t) 



+ ^ (50) 



At any time t the expectation value of the number operator for pions (in each fc-mode) 

is 

,,, ,,,, Tr<fc(t)-a.,,(t)p(0) Tr<,(0)-a^.,(0)p(t) 

<^-^^^)^ = T^^M^^ = T^) • ^''^ 

After some straightforward algebra we find 

{N^,k){t) = (2 I J-+,fc(t,to) r -l)(Ar,,fc)(0) + (I ^+,fc(t,to) r -l) . (52) 

The first term represents the "induced" and the second term the "spontaneous" particle 
production. In this approximation, particle production is a consequence of parametric am- 
plification. The Hartree-Fock states are examples of squeezed states, and the density matrix 
is a "squeezed" density matrix. The squeeze parameter (the ratio of the kernels at a time t 
to those at time t = 0) is time dependent and determines the time evolution of the states 
and density matrix. The relation between squeezed states and pion production has been 



advanced by KoganpOf although not in the context of an initial thermal density matrix. 

Thus far we have established the formalism to study the non-equilibrium evolution during 
the phase transition. A question of interpretation must be clarified before proceeding further. 
Our description, in terms of a statistical density matrix, describes an isospin invariant mixed 
state, and thus does not prefer one isospin direction over another. A real experiment will 
furnish one realization of all the available states mixed in the density matrix. However, if 
the pion correlation functions become long ranged (as a statistical average) it is clear that 
in a particular realization, at least one isospin component is becoming correlated over large 
distances, thus it is in this statistical sense that our results should be understood. 

This concludes our discussion of the formalism we will use to study the non-equilibrium 
evolution of the pion system. We now turn to a numerical analysis of the problem. 

3 Numerical Analysis 

The phenomenological set of parameters that define the linear sigma model as an effective 
low energy theory are as follows (we will be somewhat cavalier about the precise value of 
these parameters as we are interested in the more robust features of the pion correlations) 
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M„ ^ 600 MeV ; /^ ^ 95 MeV ; A ^ 4.5 

h ^ (120 MeV f ; Tc^ 200 MeV . (53) 

The above value of the critical temperature differs somewhat from the lattice estimates 
(Tg ^ 150 MeV ), and our definition of A differs by a factor four from that given elsewhere [p, 

The first thing to notice is that this is a strongly coupled theory, and unlike our previous 
studies of the dynamics of phase transitions ||12|, |13] we expect the relevant time scales to be 
much shorter than in weakly coupled theories. 

We must also notice that the linear sigma model is an effective low-energy cutoff theory. 
There are two physically important factors that limit the value of the cutoff: i) this effective 
theory neglects the influence of the nucleons, and in the Hartree approximation, the vector 
resonances are missed. These two features imply a cut-off of the order of about 2 GeV; ii) 
the second issue is that of the triviality bound. Assuming that the value of the coupling is 
determined at energies of the order of M„, its very large value implies that the cutoff should 
not be much larger than about 4-5 GeV, since otherwise the theory will be dangerously 
close to the Landau pole. From the technical standpoint this is a more important issue 
since in order to write the renormalized equations of motion we need the ratio between 
bare and renormalized couplings. In the Hartree approximation this is the "wave function 
renormalization constant" for the composite operator vf^. 

Thus we use a cutoff A = 2 GeV . The issue of the cutoff is an important one since 
(vf^) requires renormalization, and in principle we should write down renormalized equations 
of motion. In the limit when the cutoff is taken to infinity the resulting evolution should 
be insensitive to the cutoff. However the chosen cutoff is not very much larger than other 
scales in the problem and the "renormalized" equations will yield solutions that are cutoff 
sensitive. However, this sensitivity will manifest itself on distance scales of the order of 0.1 
fm or smaller, and we are interested in detecting correlations over many fermis. The size 
of the correlated regions and the time scales for their growth will be determined by the 
long wavelength unstable modes|T^ (see below), and thus should be fairly insensitive to the 
momentum scales near the cutoff. The short distance features of the correlation functions, 
such as e.g. the amplitude of the fluctuations will, however, be rather sensitive to the cutoff. 

The most severe ultraviolet divergence in the composite operator vr^ is proportional to A^. 
This divergence is usually handled by a subtraction. We will subtract this term (including 
the temperature factors) in a renormalization of the mass at t = 0. Thus 
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m|(t = 0) + 4A(7r2)(t = 0; T) = m|(t = 0; T) (54) 

where we made explicit the temperature dependence of tt^ and m^{t = 0;T). Furthermore 
we will parametrize the time dependent mass term as 



Kit) = -f 



2^2 



exp 



-2- 






M2 



2^2 



0(-t). 



(55) 



This parametrization incorporates the dynamics of the expansion and cooling processes 
in the plasma in a phenomenological way. It allows for the system to cool down with an 
effective temperature given by: 



^e//(t) =Texp 



t 



(56) 



where T is the initial value of the temperature in the central rapidity region, and t^ is a 
relaxation time. This parametrization also allows us to study a "quench" corresponding to 
the limiting case tr = 0. 

It is convenient to introduce the natural scale fm ^^ ^ 200 MeV = Mp and define the 
following dimensionless variables 



0(t) = Mpfit) ; ^^,fc(t) 



k = Mpq ; t 



T 



) f-r 



Mr) 

Tr 



Mf 



X 



Mf 



In these units 



A 

Ml 



10 ; i^ = 3 ; /7 = i^0.22 



M, 



W^ = q' + 



Ml 
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rp2 



- 1 



+ 4A/2(0) 



Thus we have to solve simultaneously the Hartree set of equations: 



(Pf , 9 



+ 7;f 



dT^ 2 



7^2 
c 



exp 



T 



-1 



+ 4A/^ + 4/AS(0,r)-if 



■ rf2 „ 9 
^ + ^+2 



T2 
^exp 



r 
-2 — 



+ 4A/^(r)+4AS(0,r)U,(r) 








(57) 

(58) 
(59) 

(60) 
(61) 



S(z,r) = (7f(f,t)-7f(0,t))/M| 

10 sin(g2;) ,. . , . ,2 



47r2 



z 



ijq{0) Hcoth 



w^„ 



10T( GeV 



(62) 
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with the boundary conditions 



m = /o ; ^ = (63) 

^l;,{0) = -^ ; ^(0) = t^, (64) 



Finally, once we find the Hartree mode functions, we can compute the total pion density as 
a function of time: 

n ( fm )3 27r2 fm ^ Jo 



dq g'(iV.,,(r)) (65) 



with (A^7r,g(T)) given by ( |52| ) in terms of the dimensionless variables. 

The mechanism of domain formation and growth is the fast time evolution of the unstable 



modes fT^ 



Let us consider first the case of if = 0. Then / = is a fixed point of the evolution 
equation for / and corresponds to cooling down from the symmetric (disordered) phase in the 
absence of explicit symmetry breaking perturbations. Let us consider the simpler situation 
of a quench (r^ = 0). The equation for the Hartree mode functions (|6TD shows that for 
g2 < 9/2 the corresponding modes are unstable at early times and grow exponentially. 

This growth feeds back on S(t) which begins to grow and tends to overcome the insta- 
bility. As the unstable fluctuations grow, only longer wavelengths remain unstable, until the 
time when 4AZI(r) ^ 9/2, at which point no wavelength is unstable. The modes will continue 
to grow however, because the derivatives will be fairly large, but since the instabilities will 
be overcome beyond this time the modes will have an oscillatory behavior. 

We expect then that the fluctuations will grow during the time for which there are 
instabilities. This time scale depends on the value of the coupling; for very small coupling, 
E(r) will have to grow for a long time before 4AS(r) ^ 9/2 and the instabilities are shut-off. 
On the other hand, for strong coupling this time scale will be rather small, and domains will 
not have much time to grow. 

It is clear that allowing for a non-zero magnetic fleld or /(O) 7^ will help to shut off the 
instabilities at earlier times, thus making the domains even smaller. 

For typical relaxation times (^1 — 2 fm/c ), domains will not grow too large either 
because the fast growth of the fluctuations will catch up with the relaxing modes and shut- 
off the instabilities (when Teff{T) < Tc) on short time scales. Thus we expect that for a 
non-zero relaxation time Tr ~ fm/c , domains will not grow too large either because the 
fluctuations will shut-off the instabilities (when T^fflr) < T^) on shorter time scales (this 
argument will be conflrmed numerically shortly). 
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We conclude from this analysis that the optimal situation for which large DCC regions 
can grow corresponds to a quench from the symmetric phase in the absence of a magnetic 
field. 

Figure (l.a) shows 4AS(r) vs r, figure (l.b) shows Ti{z,t = 1;3;5) vs z and and figure 
(l.c) shows Nt^{t) vs t for the following values of the parameters: Tc = 200 MeV ; T = 
250 MeV ; Tr = ; /(0) = 0; H = ; X = 4.5 corresponding to a quench from the 
symmetric phase at a temperature slightly above the critical temperature and no magnetic 
field. 

We clearly see that the fluctuations grow to overcome the instability in times ^ 1 f m 
and the domains never get bigger than about ~ 1.5 fm . Figure (l.c) shows that the number 
of pions per cubic fermi is about 0.15 at the initial time (equilibrium value) and grows to 
about 0.2 in times about 1-2 fermis after the quench. This pion density is thus consistent 
with having only a few pions in a pion-size correlation volume. 

Figures (2.a-c) show again the same functions but now we let the system relax from 
an initial temperature T = T^ with a relaxation time tr = I i^a jc. Notice that now the 
fluctuations grow less rapidly as they are modulated by the relaxation time, but again on a 
time scale of the order of a fm/c, they become big enough to shut-off the long wavelength 
instabilities. Figure (2.b) shows S(2;, r = 1;3;5;7) vs z. Once again, pions are correlated 
over distances of the order of a Fermi. The reason the fluctuations grow so quickly and thus 
shut off the growth of the unstable modes so quickly is the strongly coupled nature of the 
theory. 

The possibility of long range correlations exists if the initial state is in equilibrium at 
the critical temperature. In this situation there are already long range correlations in the 
initial state that will remain for some time as the temperature factors enhance the con- 
tributions for long wavelength modes since the Boltzmann factor ^ 1/k for long wave- 
length fluctuations. Figures (3.a-c) show this situation for the values of the parameters 
T = T^ = 200 MeV ; r^ = ; /(O) = ; H = ; A = 4.5. Figure (3.b) shows 
S(z,r = 1;2) vs z. In this case the number of pions per cubic fermi in the initial state is 
^ 0.12 and reaches a maximum of about 0.17 within times of the order of a fermi/c. The 
pions, however, are correlated over distances of about 4 — 5 fm with a large number of pions 
per correlation volume ~ 50. These large correlation volumes are a consequence of the initial 
long range correlations. This is the situation proposed by Gavin, Gocksch and Pisarski for 
the possibility of formation of large domains, as there is a "massless" particle in the initial 
state. 

We believe that this situation is not very likely as the central rapidity region must remain 
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in equilibrium at (or very close) to the critical temperature before the quench occurs. 

To contrast this situation with that of a weakly coupled theory, figures (4.a-c) show the 
same functions with the following values of the parameters A = 10~^ ; Tc = 200 MeV ; T = 
250 MeV ; r^ = ; /(O) = ; H = 0. Now the fluctuations are negligible up to 
times of about 4 fm/c, during which time the correlation functions grow (figure 4.b) and 
pions become correlated over distances of the order 4-5 Fermis. As shown in figure (4.c) 
the number of pions per cubic fermi becomes enormous, a consequence of a large paramet- 
ric amplification. In this extremely weakly coupled theory, the situation of a quench from 
above the critical temperature to almost zero temperature does produce a large number of 
coherent pions and domains which are much larger than typical pion sizes. This is precisely 
the situation studied previously |jl2|] within a different context. 

We have analyzed numerically many different situations in the strongly coupled case 
(A ~ 4 — 5) including the magnetic field and letting the expectation value of the sigma field 
"roll-down" etc, and in all of these cases in which the initial temperature is higher than the 
critical (between 10 — 20% higher) we find the common feature that the time and spatial 
scales of correlations are ^ 1 ira . Thus it seems that within this approach the strongly 
coupled linear sigma model is incapable of generating large domains of correlated pions. 

4 Discussions and Conclusions: 

Our study differs in many qualitative and quantitave ways from previous studies. In partic- 
ular we incorporate both quantum and thermal fluctuations and correlations in the initial 
state. In previous studies it was argued that because one is interested in long-wavelength 
fluctuations these may be taken as classical and the classical evolution equations (with cor- 
relations functions replaced by spatial averages) were studied. We think that it is important 
to quantify why and when the long-wavelength fluctuations are classical within the present 
approximation scheme. This may be seen from the temperature factors in the Hartree prop- 
agators. These are typically (incorporating now the appropriate powers of h): 



/icoth 



2kBT 



Thus, modes with wavelength k and energies Uk are classical when huji^ <^ fc^T and yield a 
contribution to the propagator 

h coth ^ 2kBT/uJk 

[2kBT\ ' 

16 



(notice the cancellation of the h). For long- wavelength components this happens when 

because the "thermal mass" (squared) for the excitations in the heat bath is -i^\T'^ /T"^ — 1]. 
For the phenomenological values of Ma- and Tc-, the "classical" limit is obtained when 

[L-l]«0.1, (66) 

c 

that is, when the initial state is in equilibrium at a temperature that is extremely close to 
the critical temperature. This is the situation that is shown in figures (3.a-c) where, indeed, 
we obtain very large correlated domains that were already present in the initial state after 
a quench from the critical temperature all the way to zero temperature. 

After an energetic collision it seems rather unlikely that the central region will be so 
close to the critical temperature. If the temperature is higher than critical, in order for 
the system to cool down to the critical temperature (or very near to it) and to remain 
in local thermodynamic equilibrium, very long relaxation times are needed, as the long- 
wavelength modes are typically critically slowed down during a transition. Long relaxation 
times will allow the fluctuations to shut off the instabilities as they begin to grow and the 
system will lose its long range correlations. This was the original argument that discarded 
an equilibrium situation as a candidate for large domains. Furthermore, typical heavy ion 
collisions or high energy processes will not allow long relaxation times (typically of a few 
fermis/c). Thus we believe that in most generic situations, a classical approximation for the 
long wavelength modes is not reliable in the Hartree approximation. We should make here a 
very important point. We are not saying that large coherent fluctuations cannot be treated 
semiclassically. They can. What we are asserting with the above analysis is that within 
the Hartree approximation, long wavelength excitations cannot be treated as classical. The 
Hartree approximation in the form used by these (and most other) authors does not capture 
correctly the physics of coherent semiclassical non-perturbative conflgurations. 

Thus although the most promising situation, within the model under investigation, is a 
quench from the critical temperature (or very close to it) down to zero temperature, it is our 
impression that this scenario is physically highly unlikely. 

There is another very tantalizing possibility for the formation of large correlated pion 
domains within the linear sigma model and that is via the creation of a critical droplet 
that will complete the phase transition (flrst order in this case) via the process of thermal 
activation over a "free energy" barrier. The small magnetic field (resulting from the small 
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up and down quark masses) introduces a small metastability [pl|] . The classical equations of 
motion allow for a solution in which the pion field is zero everywhere and a droplet in the 
sigma field (this is the 0(3) symmetric bounce responsible of thermal activation in scalar 



metastable theories in three dimensions p^] ). Using a spherically symmetric ansatz for a 
sigma droplet of radius R 

c^d(r)^/.tanh[M,(r-i?)], (67) 

and assuming, for the sake of argument, that the thin-wall approximation is reliable, we 
obtain an approximate form for the energy of the droplet 

4:71 

E ^ AnR^f^M^ - —R'hU (68) 

The critical radius is thus (this approximation is clearly reliable only as an order-of-magnitude 
estimate) 

i?c ~ 3 - 5 fm . 

By considering the fluctuations of the pions around this configuration, it is conceivable 
(although we cannot provide a more convincing argument at this stage) that the unstable 
mode of the droplet (dilation) that makes the droplet grow to complete the phase transition 
via thermal activation, produces a large amount of correlated pions. This scenario, however, 
requires supercooling (the false vacuum to be trapped) which again requires long relaxation 
times (again unlikely for strong coupling). 

As argued above, this possibility cannot be studied via a Hartree approximation which 
only provides a (select) resummation of the perturbative expansion and is probably reliable 
only for short times, before non-perturbative configurations truly develop. 

Thus we conclude that although our analysis provides a negative answer to the question 
of the possibility of large correlated domains near the chiral phase transition, these results 
are valid only within the Hartree approximation of the linear sigma model. There are several 
conceivable possibilities that would have to be studied thoroughly before any conclusions are 
reached: i) perhaps the linear sigma model is not a good candidate for studying the dynamics 
of the chiral phase transition (although it describes the universality class for the static 
properties) ii) there are large coherent field configurations (droplets) that are not captured 
in the Hartree approximation. This possibility is rather likely and is closer to the scenario 
envisaged by Bjorken, Kowalski and Taylor 0. These semiclassical coherent configurations 
may be responsible for large regions of correlated pions. An important ingredient in this 
latter case must be a deeper understanding of the dynamical (relaxation) time scales, for 
which a deeper understanding of the underlying strongly interacting theory is needed. A 



particularly relevant question is whether such a strongly coupled theory can yield enough 
supercooling so as to produce such a configuration. 

We believe that these two possibilities must be studied further to give an unequivocal 
answer to the question of large correlated domains. We are currently studying the second 
possibility in more detail. 
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Figure Captions: 
Figure l.a: 

4AS(r) vs r. T^ = 200 MeV ; T = 250 MeV ; r^ = ; /(O) = ; H = ; A = 
4.5 

Figure l.b: Ti{z, r = 1; 3; 5) vs z for the same values of the parameters as in figure (l.a) 
larger values of time correspond to larger amplitudes at the origin. 

Figure l.c: Nt^{t) vs r for the same values of the parameters as in figure (l.a). A^7r(0) = 
0.15. 

Figure 2. a: 

4AS(r) vs r. T = 250 MeV T^ = 200 MeV ; r^ = 1 ; /(O) = ; H = ; A = 4.5 

Figure 2.b: Ti{z, r = 1; 3; 5; 7) vs z for the same values of the parameters as in figure 
(2. a) larger values of time correspond to larger amplitudes at the origin. 

Figure 2.c: Nt,-{t) vs r for the same values of the parameters as in figure (2. a). A^7r(0) = 
0.15. 

Figure 3. a: 

4AS(r) vs r. T = Tc = 200 MeV ; r^ = ; /(O) = ; H = ; A = 4.5 

Figure 3.b: S(z, r = 1; 2) vs z for the same values of the parameters as in figure (3. a) 
larger values of time correspond to larger amplitudes at the origin. 

Figure S.c: NT^{r) vs r for the same values of the parameters as in figure (3. a). A^7r(0) = 
0.12. 

Figure 4. a: 

4AS(r) vs T. Te = 200 MeV ; T = 250 MeV ; r^ = ; /(O) = ; H = ; A = 
10-6 

Figure 4.b: S(z, r = 4; 6) vs z for the same values of the parameters as in figure (4. a) 
larger values of time correspond to larger amplitudes at the origin. 

Figure 4.c: Nt,-{t) vs t for the same values of the parameters as in figure (4. a). A^7r(0) = 
0.15. 
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